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Integrating  PEM  fuel  cells  effectively  with  liquid  hydrocarbon  reforming  requires  careful  system  analysis 
to  assess  trade-offs  associated  with  H2  production,  purification,  and  overall  water  balance.  To  this  end, 
a  model  of  a  PEM  fuel  cell  system  integrated  with  an  autothermal  reformer  for  liquid  hydrocarbon  fuels 
(modeled  as  C12H23)  and  with  H2  purification  in  a  water-gas-shift/membrane  reactor  is  developed  to 
do  iterative  calculations  for  mass,  species,  and  energy  balances  at  a  component  and  system  level.  The 
model  evaluates  system  efficiency  with  parasitic  loads  (from  compressors,  pumps,  and  cooling  fans), 
system  water  balance,  and  component  operating  temperatures/pressures.  Model  results  for  a  5-kW  fuel 
cell  generator  show  that  with  state-of-the-art  PEM  fuel  cell  polarization  curves,  thermal  efficiencies  >30% 
can  be  achieved  when  power  densities  are  low  enough  for  operating  voltages  >0.72  V  per  cell.  Efficiency 
can  be  increased  by  operating  the  reformer  at  steam-to-carbon  ratios  as  high  as  constraints  related  to 
stable  reactor  temperatures  allow.  Decreasing  ambient  temperature  improves  system  water  balance  and 
increases  efficiency  through  parasitic  load  reduction.  The  baseline  configuration  studied  herein  sustained 
water  balance  for  ambient  temperatures  <35  °C  at  full  power  and  <44  °C  at  half  power  with  efficiencies 
approaching  ~27  and  ~30%,  respectively. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

While  significant  work  has  been  done  on  modeling  individual 
components  of  a  PEM  fuel  cell  system,  integrated  system  sim¬ 
ulation  (including  balance-of-plant)  presents  a  challenge  for  the 
design  of  actual  fuel  cell  power  generators.  This  challenge  for  sim¬ 
ulating  balance-of-plant  is  substantially  more  complicated  when 
liquid  hydrocarbon  fuel  processing  and  subsequent  H2  purification 
must  be  incorporated  into  the  analysis.  The  needs  for  recovering 
waste-heat  and  for  sustaining  water  balance  in  both  the  fuel  pro¬ 
cessor  and  fuel  cell  stack  require  significant  thermal  integration  and 
mass  flow  recycling.  Consequently,  it  is  difficult  to  evaluate  single 
component  operation  outside  of  the  context  of  the  entire  system. 
Integrated  system  models  must  include  the  PEM  fuel  cell  stack,  the 
fuel  processing  system  (with  steam  and/or  air  input),  a  hydrogen 
purification  process,  and  the  requisite  balance-of-plant  compo¬ 
nents  (liquid  pumps,  compressors,  heat  exchangers,  and  cooling 
loops).  Detailed  component  models  with  multi-dimensional  grids, 
such  as  for  PEM  fuel  cell  stacks  [1-5]  and  for  catalytic  fuel  pro¬ 
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cessors  [6,7],  become  computationally  prohibitive  for  integrated 
system-level  analysis  for  a  broad  range  of  operating  conditions  and 
design  modifications.  As  a  result,  integrated  system  models  with 
lumped  analysis  for  individual  components  [8-12]  have  been  used 
for  evaluating  system  design  and  feasibility  of  complex  integrated 
PEM  fuel  cell  power  plants  with  hydrocarbon  fuel  processors. 

Even  simplified  component  level  analyses  for  fuel  cells  with  liq¬ 
uid  hydrocarbon  fuel  processing  show  a  wide-degree  of  variability 
in  the  detail  of  their  approaches.  Approaches  for  evaluating  sys¬ 
tem  water  balance  that  include  the  fuel  processor  and  PEM  fuel 
cell  have  been  presented  for  idealized  fuel  processors  either  with¬ 
out  consideration  of  CO  contamination  in  the  fuel  cell  or  with 
preferential  oxidation  reactors  for  CO  clean-up  [8,9].  The  previ¬ 
ous  studies  show  the  importance  of  exhaust  water  recovery  for 
maintaining  water  balance,  but  these  models  presented  in  the 
literature  have  not  fully  incorporated  the  effects  of  fuel  cell  per¬ 
formance  and  balance-of-plant  parasitic  losses  on  overall  water 
demand.  Simulating  the  integrated  performance  of  an  autothermal 
fuel  processor  with  a  PEM  fuel  cell  system  requires  proper  account¬ 
ing  of  the  thermal  integration  for  the  fuel  processor  due  to  the 
wide  variability  in  exothermicity/endothermicity  of  fuel  processors 
with  variations  in  steam  to  carbon  (S/C)  and  oxygen  to  carbon  (O/C) 
feed  ratios.  Studies  from  Argonne  National  Laboratory  have  focused 
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Nomenclature 

Acel 

geometric  area  of  single  electrochemical  cell  in  fuel 
cell  stack  (cm2) 

^H2,comb 

lower  enthalpy  of  combustion  of  H2  (J  kg-1 ) 

^HC.comb 

lower  enthalpy  of  combustion  of  hydrocarbon 

(J  kg-1 ) 

i 

average  effective  current  density  per  unit  area  in  the 
fuel  cell  stack  (A cm-2) 

^H2,mem 

mass  flow  rate  of  H2  out  of  low  pressure  side  of  the 
Pd-alloy-membrane  (kg  s-1 ) 

^HC.in 

mass  flow  rate  of  hydrocarbon  into  fuel  reformer 
(kgs-1) 

^cells 

number  of  electrochemical  cells  in  fuel  cell  stack 

O/C 

oxygen  (from  02)  to  carbon  (from  fuel)  atomic  ratio 
in  fuel  reformer  inlet 

S/C 

steam  to  carbon  (from  fuel)  ratio  in  fuel  reformer 
inlet 

^amb 

ambient  temperature  (°C  or  K) 

^exh 

exhaust  condenser  outlet  temperature  (°C  or  K) 

V'cell 

voltage  across  an  individual  electrochemical  cell  in 
stack  (V) 

Wlost 

parasitic  power  lost  due  to  balance-of-plant  compo¬ 
nents  and  dc-ac  inverter  (W) 

Wnet 

net  electrical  power  out  (W) 

Greek  Symbols 

*7bop 

effective  efficiency  associated  with  power  delivered 
to  power  produced  by  stack 

??fc 

effective  efficiency  of  fuel  cell  stack 

Vfp 

effective  efficiency  of  fuel  processor  and  Pd-alloy- 
membrane  purifier  combined 

Vth 

overall  effective  system  efficiency 

on  liquid  hydrocarbon  reforming  without  H2  purification  [8]  and 
with  a  preferential  oxidation  reactor  for  CO  clean-up  and  an  ide¬ 
alized  reformate  tolerant  stack  [10,12].  Other  studies  investigated 
a  detailed  model  of  a  down-the-channel  reformer  with  a  Pd-based 
H2  separation  membrane  [11  ],  similar  to  what  is  considered  in  the 
present  paper.  Nonetheless,  previous  studies  on  liquid  hydrocarbon 
processing  did  not  fully  address  the  impact  of  balance-of-plant  on 
fuel  processor  loads,  fuel  cell  power  demand,  and  system  temper¬ 
atures,  particularly  for  small  portable  power  (<10  kW)  applications 
which  underlie  the  current  study. 

The  current  study  attempts  to  extend  the  previously  published 
analyses  of  integrated  PEM  fuel  cell  systems  with  liquid  hydro¬ 
carbon  fuel  processors  by  performing  mass  and  energy  balances 
on  a  workable  system  that  includes  an  integrated  water-gas-shift 
reactor/H2  palladium  (Pd)  alloy  membrane  separator  to  ensure 
reliable  PEM  performance.  The  present  analysis  is  for  an  inte¬ 
grated  system  that  produces  between  0.5  and  5  kW  net  power  out 
using  a  kerosene  surrogate  (Ci2H23  with  thermodynamic  proper¬ 
ties  defined  previously  [13]).  By  exploring  the  mass  and  energy 
balances  of  the  fuel  processor,  H2  purification,  fuel  cell  system,  and 
balance-of-plant,  this  model  provides  critical  assessment  not  only 
of  the  important  trade-offs  in  terms  of  thermal  efficiency  and  water 
balance,  but  also  of  the  necessary  component  performance  require¬ 
ments  and  range  of  operating  conditions  amenable  for  efficient  and 
sustainable  system  operation.  The  following  sections  will  present 
the  system  layout  and  modeling  approach,  followed  by  a  discussion 
of  simulation  results  on  critical  system  performance  trends  for  a 
range  of  operating  conditions  as  well  as  component  performance 
criteria. 


2.  System  model 

2.1.  System  configuration 

A  liquid-fuel  processor  and  H2  purification  system  can  be  inte¬ 
grated  into  a  PEM  fuel  cell  power  generator  in  a  number  of  possible 
flow  configurations,  each  with  different  approaches  to  thermal  and 
water  management.  System-level  modeling,  as  presented  in  this 
study  and  earlier  work  [8,10,11],  is  needed  to  compare  an  array  of 
possible  configurations.  Furthermore,  assessing  the  viability  and 
robustness  of  a  particular  design  requires  the  determination  of  per¬ 
formance  over  a  range  of  working  conditions.  The  paper  focuses  on 
one  particular  configuration,  illustrated  in  Fig.  1,  that  integrates  a 
PEM  fuel  cell  stack  with  a  liquid  fuel  processor,  a  water-gas-shift/Pd 
membrane  reactor  for  H2  purification,  and  additional  balance-of- 
plant  subsystems. 

The  H2  supply  for  the  fuel  cell  anode  is  obtained  from  an 
autothermal  fuel  reformer,  followed  by  a  water-gas-shift  reactor 
that  is  integrated  with  a  Pd-alloy-membrane  for  H2  separation. 
Because  the  fuel  processor  is  upstream  of  the  H2  separation  mem¬ 
brane,  it  must  operate  at  sufficient  pressure  to  provide  a  pressure 
differential  (in  this  case,  6.0barg)  for  driving  H2  across  the  mem¬ 
brane.  The  model  assumes  a  desulfurized  fuel  source;  and  thus  no 
pressure  drop  penalty  for  sulfur  removal  is  incorporated  in  the  fuel 
processor.  Pumps  for  fuel  and  water,  along  with  a  high-pressure, 
low-flow-rate  compressor  for  air,  provide  inlet  flows  to  the  fuel 
processor.  On  the  low-pressure  side  of  the  Pd-alloy-membrane, 
purified  H2  is  carried  from  the  membrane  by  a  flow  of  sweep 
steam.  The  sweep  steam  also  acts  as  a  diluent,  reducing  H2  par¬ 
tial  pressure  and  promoting  diffusion  across  the  membrane.  The 
fraction  of  H2  not  diffusing  across  the  membrane,  along  with  all 
other  reformate  species,  exits  the  membrane  separator  as  reten- 
tate. 

The  H2/H20  mixture  flow  exiting  the  low-pressure  side  of  the 
membrane  separator  cools  as  it  mixes  with  H2 /water  vapor  recir¬ 
culated  from  the  anode  exhaust,  and  this  combined  flow  is  further 
cooled  through  a  heat  exchanger  before  entering  the  anode.  Part  of 
the  water  vapor  condenses  in  the  heat  exchanger  and  is  available 
for  re-use  in  the  system.  Additional  liquid  water  is  recovered  from 
flow  exiting  the  anode.  Collected  water  that  is  needed  for  sweep 
steam  is  pumped  to  a  steam  generator  and  vaporized.  Any  surplus 
liquid  water  can  be  stored  or  sent  to  the  fuel  processor  as  required. 
Flow  through  the  fuel  cell  anode  is  recirculated,  ensuring  uniform 
H2  distribution  and  preventing  liquid  water  accumulation  along  the 
anode  channels  while  utilizing  all  available  H2.  A  compressor  pro¬ 
vides  the  pressure  differential  needed  to  induce  flow  recirculation 
through  the  anode.  N2  buildup  in  an  anode  recirculation  due  to  leak¬ 
age  across  the  Nation  membrane  will  require  some  purging,  and  a 
recent  study  indicates  that  purging  to  maintain  N2  content  below 
1.5%  costs  the  fuel  cell  stack  efficiency  approximately  1%  [14].  Such 
a  purge  cycle  is  not  readily  modeled  in  a  steady-state  model  and 
the  small  impact  on  overall  efficiency  in  the  system  motivates  the 
assumption  in  this  study  that  the  purge  in  the  anode  recirculation 
loop  is  negligible. 

Air  is  supplied  to  the  fuel  cell  cathode  by  a  rotary  vane  com¬ 
pressor  and  is  cooled  through  a  heat  exchanger  after  exiting  the 
compressor.  A  gas-to-gas  humidifier  with  a  Nafion-based  water 
membrane  humidifies  the  air  stream  entering  the  cathode  with 
water  removed  from  the  cathode  exhaust. 

Both  retentate  and  cathode  exhaust  are  fed  to  a  catalytic  burner. 
Burner  exhaust  provides  heat  input  to  the  reformer  pre-heater, 
reformer  steam  generator,  and  reformer  fuel  vaporizer.  Further 
downstream,  the  exhaust  is  cooled  in  a  condenser  to  recover  water. 
The  water  thus  recovered  can  be  recycled  to  the  fuel  processor,  but 
cannot  be  used  for  sweep  steam  or  cathode  humidification  because 


1058 


J.B.  Pearlman  et  al.  /  Journal  of  Power  Sources  185  (2008)  1056-1065 


_3 

_4 

_5 

_7 

10 

= 

12 

13 

14 

15 


Fuel  supply  to  burner 

Autothermal  reformer  to  sweep  steam  generator  (ref  exh  flow) 
Sweep  steam  generator  to  water  gas  shift  reactor  (ref  exh  flow) 
Retentate  to  exhaust  burner 

Exhaust  burner  to  reformer  pre-heater  (burner  exhaust  flow) 
Reformer  pre-heater  to  reformer  steam  generator  (burner  exh  flow) 
Reformer  steam  generator  to  fuel  vaporizer  (burner  exh  flow) 

Fuel  vaporizer  to  condenser  (burner  exh  flow) 

Condensed  water  to  reformer  water  reservoir 
Condensate  from  anode  recirculation  loop  to  sweep  water  reservoir 
Coolant  radiator  to  LP  compressor  heat  exchanger  (cooling  loop) 
Coolant  radiator  to  anode  inlet  heat  exchanger  (cooling  loop) 
Coolant  radiator  to  stack  cooling  (cooling  loop) 

LP  compressor  heat  exchanger  to  coolant  pump  (cooling  loop) 
Anode  inlet  heat  exchanger  to  coolant  pump  (cooling  loop) 


Liquid  fuel 

Air  supply  and  cathode  stream 

Vaporized  fuel  /  water  /  air  mixture 

Reformate  and  retentate 

Liquid  water  and  steam  (fuel  processor) 

Liquid  water  and  steam 

Permeate  and  anode  stream 

Liquid  coolant 

Products  from  burner 


Fig.  1.  Schematic  diagram  showing  realization  of  PEM  fuel  cell  system  integrated  with  liquid  fueled  ATR  reactor,  water-gas-shift/Pd-alloy-membrane  reactor,  and  retentate 
exhaust  burner  for  waste-heat  recovery. 


of  the  presence  of  H2C03  and  other  potential  impurities  that  may 
impact  long-term  stability  of  the  PEM  fuel  cell  stack. 

A  liquid  coolant  loop  removes  heat  from  the  fuel  cell  stack,  in 
addition  to  cooling  heat  exchangers  for  the  anode  recirculation  flow 
and  cathode  inlet  flow.  The  coolant  is  circulated  through  a  radiator, 
where  the  accumulated  heat  is  rejected  to  the  ambient. 

2.2.  System  modeling  approach 

The  model  data  is  imported  into  MS  Excel  with  a  complex  Visual 
Basic  program  running  as  a  macro  that  determines  steady-state 
operating  conditions  and  system  performance  through  a  nested 
iterative  sequence.  The  recycling  of  mass  and  heat  flows,  along  with 
the  feedback  between  system  operating  conditions  and  parasitic 
loads,  requires  such  an  iterative  solution  technique.  A  given  iter¬ 
ation  works  by  using  current  values  for  flow  rates  and  parasitic 
power  losses  for  the  system  to  determine  a  new  value  for  power 
produced  by  the  fuel  cell  stack.  This  power  demand  is  then  used 


to  recalculate  fuel  and  air  flows,  thermal  loads  for  cooling  compo¬ 
nents,  and  associated  parasitic  loads  (from  compressors,  pumps, 
and  cooling  fans)  for  system  operation.  Fuel  flow  rate  to  the  sys¬ 
tem  is  determined  iteratively  through  a  species  balance  to  satisfy 
H2  demand  from  the  Pd-alloy-membrane  purifier.  The  new  flow 
rates  thus  are  used  to  update  the  parasitic  loads,  and  gross  power 
demand  for  the  fuel  cell  stack  is  calculated  from  the  sum  of  net  (i.e., 
usable)  power  demand  and  the  parasitic  loads. 

An  iterative  solution  technique  is  used  to  determine  tempera¬ 
tures  in  each  component  to  satisfy  entropy  constraints  and  energy 
balances.  Successively  repeated  iterations  that  solve  species  and 
energy  balance  equations  are  conducted  until  the  following  criteria 
were  met.  ( 1 )  The  difference  in  overall  gross  power  demand  and  fuel 
cell  power  produced  are  within  10-4  W.  (2)  The  internal  enthalpy 
balance  for  each  component  is  within  10-3  W.  (3 )  The  molar  balance 
for  the  fuel  cell  anode  recirculation  loop  is  within  10-9gmols-1. 
The  iterative  solver  thus  provides  a  steady-state  solution  of  flow 
conditions  and  also  necessary  parasitic  loads  from  compressors 
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and  pumps  and  fans  for  heat  removal.  This  solution  approach  is 
similar  in  nature  to  the  approach  for  a  system-level  model  of  a 
methane-fueled  PEM  fuel  cell  system  studied  earlier  [10],  although 
that  system  relied  on  CO  clean-up  in  a  preferential  oxidation  reactor 
rather  than  a  Pd-alloy-membrane  for  H2  purification. 

Thermodynamic  properties  are  calculated  based  on  JANAF 
tables.  Thermodynamics  of  a  surrogate  kerosene  (C12H23)  is  used 
for  the  liquid  fuel  [15].  Liquid  enthalpies  are  calculated  by  subtract¬ 
ing  the  heat  of  vaporization  and  subsequent  sensible  heat  removal 
(based  on  a  constant  liquid  Cp)  from  the  gas-phase  value. 

2.3.  Component  model  description 

Most  components  as  shown  in  Fig.  1  are  simulated  with  zero- 
order  thermodynamic  models.  Exceptions  include  the  combined 
water-gas-shift/Pd-alloy-membrane  reactor,  which  is  simulated  by 
a  ID  flow  model  in  order  to  assess  the  shift  in  H2  driving  force 
across  the  membrane  as  a  function  of  length  through  the  reac¬ 
tor.  Models  for  critical  heat  exchangers  (the  exhaust  condenser, 
coolant  radiator,  anode  flow  condenser,  and  air  compressor  outlet 
heat  exchanger,  all  indicated  in  Fig.  1 )  utilize  component  geometry 
to  determine  power  needs  for  fans  and  coolant  pump  associated 
with  those  heat  exchange  processes. 

2.3.1  Fuel  cell  stack  model 

The  fuel  cell  is  modeled  as  an  isothermal  reactor  with  cell 
voltage,  Vcen,  versus  current  density,  i  (A cm-2),  defined  by  a  polar¬ 
ization  curve  comparable  to  recently  published  data  for  a  system 
based  on  Ballard  Power  Systems  technology  [13].  The  open  circuit 
voltage  is  shifted  up  or  down  based  on  changes  in  the  Nernst  poten¬ 
tial  with  changes  in  reactant  partial  pressures  and  temperatures  Tfc. 
Since  the  power  demands  for  all  cases  examined  in  this  study  drive 
the  fuel  cell  beyond  the  polarization  activation  region  of  the  Vcell-i 
curve,  a  linear  approximation  of  the  Vcell-i  curve  satisfactorily  rep¬ 
resents  the  region  of  the  polarization  curve  that  defines  the  range  of 
operation  for  the  power  conditions  explored  in  this  study.  Isother¬ 
mal  operation  of  the  fuel  cell  is  assumed  because  the  system  model 
incorporates  a  coolant  loop  with  variable  flow  rates  controlled  to 
maintain  a  near  constant  Tfc  =  65°C,  which  is  a  reasonable  value 
for  long-term  durable  operation  of  Nafion-based  membrane  elec¬ 
trode  assemblies  (MEA’s).  An  energy  balance  is  used  to  calculate 
the  coolant  flow  load  and  the  consequential  pump  and  fan  work 
required  for  the  fuel  cell  radiator  loop.  The  gas-to-gas  humidifier, 
which  recaptures  humidity  for  the  cathode  inlet  from  the  cathode 
exhaust,  is  assumed  to  be  large  enough  to  humidify  the  cathode 
inlet  flow  fully  at  its  exhaust  temperature  from  the  humidifier. 

For  the  fuel  cell  stack  model,  the  water  balance  requires  knowl¬ 
edge  of  where  the  water  product  is  removed  (either  with  the  anode 
or  cathode  flows).  Recent  detailed  two-phase  flow  MEA  models  in 
agreement  with  experiments  have  shown  that  the  net  fraction  of 
product  water  which  diffuses  back  across  the  Nation  membrane  to 
the  anode  depends  strongly  on  the  membrane  thickness,  the  H2 
stoichiometry  in  the  anode,  and  the  humidity  ratio  in  the  anode 
[5].  This  net  fraction  is  typically  quite  small  or  near  zero,  but  for 
thin  Nation  membranes  characteristics  of  the  stack  modeled  in  this 
study,  a  slight  diffusion  toward  the  anode  is  expected  if  the  anode 
is  not  fully  humidified.  In  this  study  a  net  fraction  of  15%  product 
water  back  to  the  anode  is  used  as  a  constant  stack  average  value. 
This  value  has  a  minor  impact  on  system  efficiency  and  the  net 
system  water  balance. 

2.3.2.  Fuel  processor  and  Ft2  membrane  reactor  models 

Fuel  conversion  in  the  reformer  takes  place  through  reactions 
with  both  steam  and  02.  Steam-to-carbon  (S/C)  and  oxygen-to- 
carbon  (O/C)  ratios  are  user-specified  inputs  and  set  high  enough 


for  sufficient  H20  and  02  inputs  to  achieve  complete  fuel  conver¬ 
sion.  The  exit  reformate  composition  is  based  on  water-gas-shift 
equilibrium  of  H2,  CO,  C02  and  H20  at  the  calculated  outlet  tem¬ 
perature.  This  approach  does  not  capture  the  potential  for  CH4  in 
the  reformate  stream.  CH4  concentrations  in  ATR  reformate  have 
been  shown  to  be  below  0.5%  for  operating  pressures  up  to  3  bar 
and  outlet  temperatures  of  1000  K  or  higher  at  the  S/C  ratios  in  this 
study  [  10].  A  complimentary  program  with  this  modeling  study  has 
been  developing  an  ATR  following  an  earlier  reported  design  [  16]  for 
operating  pressures  as  modeled  here,  and  results  have  shown  simi¬ 
larly  low  CH4  outlet  concentrations.  Such  low  concentrations  of  CH4 
will  not  significantly  impact  the  partial  pressure  of  H2  through  the 
membrane  reactor  nor  impede  H2  diffusion  through  the  Pd-alloy- 
membrane,  and  thus,  the  CH4  concentrations  in  the  reformate  are 
assumed  negligible  and  all  carbon  in  the  reformate  is  assumed  to 
be  in  CO  or  C02. 

For  the  cases  presented  here,  the  fuel  processor  inlet  tem¬ 
perature  Trefin  is  not  allowed  to  increase  above  500  °C,  assuming 
sufficient  heat  input  for  water  and  fuel  vaporization,  with  subse¬ 
quent  gas  flow  preheating.  While  Tref,in  is  adequate  for  all  reactor 
conditions  with  an  exothermic  reactant  mixture,  it  is  too  low  for 
conditions  where  endothermic  steam  reforming  is  the  dominant 
reaction.  The  current  study  includes  an  analysis  of  the  performance 
trade-offs  between  varying  S/C  and  O/C  ratios,  with  the  sum  of  these 
held  constant  at  2.4  and  a  mildly  exothermic  base  case  of  S/C  =  1.6 
and  O/C  =  0.8. 

The  water-gas-shift/Pd-alloy-membrane  reactor  is  modeled 
with  the  reformate  stream  and  the  sweep  stream  in  counter-flow. 
The  integrated  water-gas-shift  reaction  increases  H2  driving  poten¬ 
tial  along  the  length  of  the  Pd-alloy-membrane.  Furthermore,  as 
gases  flow  along  the  length  of  the  reactor,  H2  diffusion  across  the 
membrane  favors  additional  H2  production  in  the  reformate  stream 
to  maintain  equilibrium.  To  simulate  this  effect  more  accurately, 
the  water-gas-shift  reactor  is  discretized  into  four  segments  with 
local  energy  balances  solved  iteratively  to  determine  temperature 
and  mole  fraction  profiles  along  the  length  of  the  reactor.  In  the 
discrete  approximation,  the  water-gas-shift  reaction  is  assumed  to 
reach  equilibrium  at  the  entrance  to  each  segment. 

The  reactor  is  assumed  to  be  of  sufficient  size  that  H2  diffusion 
is  not  limited  by  mass  transfer.  Therefore,  the  amount  of  hydrogen 
that  diffuses  through  the  membrane  is  based  upon  specified  H2 
recovery  fraction.  This  represents  the  ratio  of  H2  recovered  through 
the  membrane  to  the  sum  of  H2  entering  in  the  reformate  stream 
and  H2  produced  through  water-gas-shift  in  the  membrane  reac¬ 
tor.  In  the  present  study,  membrane  H2  recovery  fraction  is  varied 
with  a  baseline  case  set  at  75%.  Significantly  higher  H2  recovery 
efficiencies  are  not  considered  feasible,  because  of  the  loss  of  H2 
partial  pressure  difference  that  serves  as  a  driving  potential  across 
the  membrane. 


2.3.3.  Compressor  and  pump  models 

The  low-pressure-air  compressor  supplies  air  to  the  fuel  cell 
cathode.  Inlet  conditions  are  set  to  ambient,  and  the  outlet  pres¬ 
sure  is  determined  to  match  the  cumulative  pressure  drop  through 
components  in  the  cathode  loop.  The  outlet  enthalpy  and  tem¬ 
perature  are  determined  from  an  isentropic  efficiency,  which  is 
determined  through  cubic  polynomials  fitted  to  publicly  available 
data  for  a  rotary  vane  compressor  [17].  These  account  for  varia¬ 
tions  in  pressure  ratio  and  mass  flow  rate.  The  required  work  input 
is  then  determined  by  the  enthalpy  change  through  the  compres¬ 
sor,  plus  any  heat  loss  to  the  ambient.  A  similar  approach  is  adopted 
for  the  high-pressure  air  compressor  and  the  anode  recirculation 
compressor  models,  except  that  isentropic  efficiency  depends  only 
on  pressure  ratio. 
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Power  required  by  the  coolant  pump,  reformer  water  pump,  fuel 
pump,  and  sweep  steam  water  pump  are  determined  from  the  isen- 
tropic  value  divided  by  a  user-specified  isentropic  efficiency.  For  the 
relatively  small  power  demand  of  the  pumps  compared  to  other 
parasitic  loads  in  the  system,  isentropic  efficiencies  are  simply  set 
constant  at  0.50. 

2.3.4.  Heat  exchanger,  fan,  and  exhaust  burner  models 

The  system  model  includes  individual  heat  exchanger  models 
for  the  exhaust  condenser,  low-pressure-air  compressor  cooler, 
anode  flow  condenser,  and  coolant  radiator.  These  are  based 
on  geometry,  with  corresponding  heat  transfer  and  pressure 
drop  correlations,  intended  for  maximum  heat  rejection  in  min¬ 
imal  volumes.  The  burner  exhaust  condenser  and  coolant  loop 
radiator  are  modeled  as  air-cooled  cross-flow  heat  exchang¬ 
ers,  while  the  air  compressor  outlet  heat  exchanger  and  the 
anode  flow  condenser  are  modeled  as  counterflow  plate  heat 
exchangers.  A  liquid  coolant  loop  provides  cooling,  via  paral¬ 
lel  branches,  to  each  plate  heat  exchanger  and  to  the  fuel  cell 
stack. 

For  the  exhaust  condenser,  inlet  temperature  and  species  mole 
fraction  values  for  the  burner  exhaust  gas  are  received  from  the 
fuel  vaporizer.  Inlet  conditions  for  the  cooling  air  are  initially  set  at 
ambient  conditions.  The  inlet  air  pressure  is  subsequently  updated 
after  the  air-side  pressure  drop  across  the  heat  exchanger  is  cal¬ 
culated.  A  target  outlet  temperature  for  the  system  exhaust  flow 
is  defined  as  10  °C  higher  than  ambient  temperature,  but  not  less 
than  30  °C  because  of  the  diminishing  benefit  in  water  recovery 
below  this  point.  The  exhaust  gas  outlet  temperature  is  normally 
held  fixed  for  given  ambient  conditions,  but  may  be  relaxed  as 
described  below.  With  the  burner  exhaust  gas  inlet  and  outlet  con¬ 
ditions  specified,  the  required  heat  load  for  the  exhaust  condenser 
is  defined  through  an  enthalpy  balance.  The  average  outlet  temper¬ 
ature  of  the  cooling  air  is  then  adjusted  iteratively  until  the  heat 
removal  capacity  of  the  heat  exchanger  matches  the  required  heat 
load.  As  the  air  outlet  temperature  is  adjusted,  it  varies  the  heat 
removal  capacity  of  the  exhaust  condenser  by  changing  the  log 
mean  temperature  difference  across  the  heat  exchanger.  Adjusting 
the  air  outlet  temperature  also  determines  the  required  air  mass 
flow  rate  to  satisfy  an  enthalpy  balance  between  burner  exhaust 
gas  and  the  cooling  air.  This  affects  the  friction  factor  and,  in  the 
case  of  turbulent  conditions,  the  heat  transfer  coefficient  through 
changes  in  flow  velocity  and  Reynolds  number.  To  prevent  second- 
law  violations,  the  solver  is  also  constrained  to  keep  the  average 
air  outlet  temperature  less  than  the  burner  exhaust  inlet  tempera¬ 
ture.  Similarly,  the  burner  exhaust  outlet  temperature  must  remain 
above  the  inlet  temperature  of  the  cooling  air. 

To  prevent  uncontrolled  growth  in  condenser  fan  power  under 
high  heat  load  conditions,  velocity  of  the  cooling  air  is  constrained 
to  a  maximum  of  3  m  s-1 .  When  this  constraint  is  active,  the  model 
relaxes  burner  exhaust  gas  outlet  temperature  to  the  extent  needed 
to  match  heat  removal  from  the  gas  with  available  heat  exchanger 
capacity. 

Constant  wall  temperature  conditions  are  assumed  for  the  pur¬ 
pose  of  determining  the  convection  coefficient  under  laminar 
conditions.  For  turbulent  flow  conditions,  the  friction  factor  inside 
the  tubes  is  determined  based  on  the  relationship  by  Petukhov 
[18]  and  the  internal  tube  heat  transfer  coefficient  is  determined 
using  Gnielinski’s  correlation  [19]  for  Nusselt  number.  In  the  transi¬ 
tion  region  between  laminar  and  turbulent  conditions,  the  friction 
factor  and  convection  coefficient  are  modeled  as  varying  linearly 
between  their  high-end  laminar  and  low-end  turbulent  values. 

The  pressure  drop  of  the  cooling  air  flow  is  determined  from  fric¬ 
tion  factor  based  on  the  correlation  of  Chang  et  al.  [20],  which  was 
found  to  predict  reasonable  values  of  friction  factor  over  the  range 


of  conditions  modeled.  For  the  air-side  heat  transfer  coefficient,  the 
model  determines  the  Nusselt  number  from  a  Chilton-Colburn  fac¬ 
tor  for  a  radiator  with  louvered  fins  using  the  correlation  found  in 
Chang  et  al.  [21  ]  for  a  Reynolds  number  (ReLp)  between  the  limits  of 
100  and  3000  which  falls  within  the  range  of  conditions  modeled 
in  this  study.  For  ReLp  <  100,  the  model  calculates  the  convection 
coefficient  based  on  a  constant  Nusselt  number,  as  evaluated  at  the 
lower  ReLP  limit. 

The  required  fan  power  is  determined  by  the  enthalpy  difference 
between  ambient  air  entering  the  fan  and  air  exiting  the  fan  at  a 
sufficient  pressure  to  match  the  pressure  drop  across  the  exhaust 
condenser.  The  exit  enthalpy  is  calculated  based  on  an  isentropic 
enthalpy  increase  divided  by  isentropic  fan  efficiency.  For  given 
combinations  of  flow  rate  and  pressure  increase,  the  isentropic  effi¬ 
ciency  is  determined  through  a  map  based  on  characteristic  curves 
measured  for  an  axial  fan.  Flow  rates  and  pressure  ratios  in  the 
map  are  expressed  relative  to  their  maximum  values,  such  that 
the  efficiency  map  can  be  scaled  to  varying  fan  sizes.  Because  of 
accuracy  limitations  of  the  efficiency  map  at  air  velocities  less  than 
0.12  ms-1,  fan  power  is  calculated  as  scaling  with  the  cube  of  air 
velocity  under  these  conditions. 

The  coolant  loop  radiator  receives  the  coolant  inlet  temperature 
and  species  mole  fraction  values  from  the  coolant  pump.  The  radi¬ 
ator  fan  speed  remains  bounded  under  all  conditions  tested  with 
no  active  constraint  on  maximum  air  velocity,  and  therefore,  no 
provision  for  relaxing  the  coolant  outlet  temperature  is  needed. 
Otherwise,  calculations  in  the  coolant  loop  radiator  are  handled 
similarly  to  those  in  the  exhaust  condenser. 

The  anode  flow  condenser  transfers  heat  from  hot  anode  gas 
to  the  coolant  loop.  The  friction  factor  and  the  turbulent-range 
heat  transfer  coefficient  for  both  flows  are  calculated  as  they  are 
for  the  burner  exhaust  flow  in  the  exhaust  condenser.  For  laminar 
conditions,  the  heat  transfer  coefficient  for  both  flows  is  based  on 
constant  surface  temperature  Nusselt  numbers  determined  by  Kays 
and  Crawford  [22]  over  a  range  of  channel  aspect  ratios.  The  coolant 
outlet  temperature  is  limited  to  a  maximum  of  90  °C  to  prevent  boil¬ 
ing.  This  has  the  effect  of  enforcing  a  minimum  coolant  mass  flow 
rate  through  the  heat  exchanger. 

The  anode  flow  condenser  is  sized  for  operation  at  design  con¬ 
ditions  that  represent  the  most  demanding  of  the  anticipated 
operating  conditions.  Under  less  severe  conditions,  the  required 
heat  load  is  smaller,  which  results  in  a  lower  mass  flow  rate 
and  greater  temperature  rise  of  the  coolant.  When  the  heat  load 
decreases  to  the  point  that  the  coolant  reaches  its  maximum  per¬ 
mitted  outlet  temperature,  the  coolant  mass  flow  rate  cannot  be 
reduced  any  more.  If  the  heat  load  is  decreased  further,  the  anode 
flow  outlet  temperature  must  drop  below  its  set  target  value  in 
order  for  enthalpy  change  in  the  anode  flow  to  match  the  heat 
removal  capacity  of  the  anode  flow  condenser.  The  model  detects 
that  this  condition  has  occurred  when  the  energy  balance  resid¬ 
ual  for  calculating  the  coolant  outlet  temperature  iterations  fails  to 
achieve  convergence  with  all  imposed  temperature  constraints.  In 
this  case,  the  model  allows  the  anode  flow  outlet  temperature  to  be 
adjusted  as  a  design  variable  and  repeats  the  iterative  calculation 
of  coolant  outlet  temperature. 

If  the  required  anode  flow  enthalpy  change  (i.e.,  heat  load)  con¬ 
tinues  to  drop,  the  anode  flow  outlet  temperature  is  allowed  to 
decrease  until  it  reaches  within  0.01  °C  of  the  coolant  inlet  temper¬ 
ature.  If  the  required  heat  transfer  to  the  coolant  decreases  beyond 
this  point,  it  becomes  insensitive  to  the  heat  exchanger  channel 
length.  Under  such  a  condition,  the  required  enthalpy  change  of 
the  anode  flow  can  take  place  in  less  than  fixed  length  of  the  heat 
exchanger  channels.  The  model  responds  by  relaxing  the  chan¬ 
nel  length  and  re-calculating  the  heat  exchanger  enthalpy  balance. 
The  channel  length  at  which  the  solver  converges  represents  the 


Table  1 

Critical  baseline  system  parameters 
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Parameter  Value 


Ambient  conditions 

Temperature,  Tamb  (I<)  303.15 

Relative  humidity,  0amb  50% 

Fuel  cell  conditions  and  properties 

Operating  temperature,  Tcen  (K)  338.15 

Pressure  drop  across  stack  at  1 A  cm-2  current  density  (bar)  0.5 

Stoichiometric  ratio  of  cathode  air  flow  1 .8 

Stoichiometric  ratio  of  anode  H2  flow  1.5 

Fraction  of  water  produced  in  fuel  cell  transported  to  anode  15% 

Number  of  cells  per  stack  25 

Stack  membrane  area  per  cell  (cm2 )  487 

Current  density  at  8000  W  gross  power,  iref  (A  cm-2 )  1 .0 

Fuel  reformer  conditions 

Inlet  oxygen  to  carbon  (O/C)  ratio  0.8 

Inlet  steam  to  carbon  (S/C)  ratio  1.6 

Maximum  inlet  temperature,  Trefjin  max  (K)  773 

Minimum  pre-heater  approach  temperature  (K)  20 

Minimum  approach  temperature  for  steam  generator  (K)  40 

Maximum  temperature  out  of  steam  generator  (K)  470 

Water-gas-shift  membrane  reactor  conditions: 

Differential  pressure  across  membrane  (bar)  6.0 

Sweep  steam  to  H2  ratio  in  the  membrane  separator  0.20 

Membrane  reactor  H2  recovery  fraction  75% 

Temperature  out  of  the  sweep  steam  generator  (K)  450 

Balance-of-plant  conditions: 

Radiator  coolant  outlet  temperature  (K)  333.15 

Burner  exhaust  condenser  temperature  difference  (K)  10 

Electric  inverter  efficiency  93% 

Electric  motor  efficiency  90% 


fraction  of  the  channel  length  actually  needed  to  accomplish  the 
required  enthalpy  change. 

No  condensation  takes  place  in  the  air  compressor  outlet  heat 
exchanger,  but  it  is  otherwise  similar  to  the  anode  flow  condenser. 
Both  heat  exchangers  use  the  same  modeling  approach. 

Exhaust  gas  from  the  retentate  burner  provides  heat  inputs  to 
the  reformer  steam  generator,  fuel  vaporizer,  and  reformer  pre¬ 
heater.  Each  of  these  is  represented  by  a  zero-order  model  that 
determines  the  burner  exhaust  outlet  temperature  iteratively  to 
satisfy  an  enthalpy  balance.  Maximum  allowable  outlet  tempera¬ 
tures  are  specified  for  the  flows  being  heated.  A  zero-order  model 
also  represents  heat  rejection  from  the  fuel  cell  stack  to  the  coolant 
loop.  In  this  model,  the  coolant  outlet  temperature  is  limited  to 
80 °C  in  order  to  prevent  excessive  stack  temperature  and  the 
accompanying  risk  of  Nation  membrane  dryout. 

3.  Results 

3.1  Baseline  system  performance 
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Fig.  2.  Overall  system  efficiency  and  total  water  balance  (in  gs-1  accumulated)  as 
functions  of  net  electric  power  for  the  liquid-fueled  PEM  fuel  cell  system  at  baseline 
conditions  (see  Table  1)  at  three  different  ambient  temperatures.  Numbers  on  top 
of  bars  indicate  Tamb. 


of-plant  ?7bop*  These  three  sub-system  efficiencies  are  shown  in 
each  expression  in  parenthesis  in  the  right-hand-side  of  the  fol¬ 
lowing  equation: 


77  th  =  77fpT7fcT7bop  = 


777H2,mem  ^H2,comb  A  l  nCe\\s  A:ell  7^cell 


77?HC,in  ^HC.comb  J  y  771h2,  mem  ^H2,  comb 


f  Ticells  ^cell  7^cell  ^lost 

y  T^cells  Acen  iVce\\ 


(2) 


?7fp  is  proportional  to  the  membrane  recovery  fraction  which  in 
this  model  is  set  constant,  and  variations  in  p^p  are  primarily  due 
to  changes  in  the  S/C  and  O/C  ratios  of  the  fuel  processor.  p¥C  is  pro¬ 
portional  to  Vcell  because  for  the  recirculated  anode  flow  at  steady 
state,  i  will  be  proportional  to  ihu2,mem-  Thus,  as  Vce\\  drops  with 
higher  i  and  higher  stack  power  densities,  p?c  will  tend  to  fall  off, 
and  this  decrease  in  with  increased  power  can  more  than  off¬ 
set  improved  efficiencies  in  the  compressors  as  they  approach  their 
capacity  with  increased  system  power.  This  trend  is  illustrated  in 
Table  2  which  shows  key  system  performance  metrics  for  baseline 
conditions  at  half-  and  full-load.  Fig.  2  shows  pth  as  a  function  of 
VVnet  for  3  different  Tam b:  10,  30,  and  50 °C.  For  all  three  Tamb,  pth 
decreases  as  Wnet  rises  above  1500  W.  This  is  due  to  the  drop  in 
Pyc-  The  decrease  in  pth  at  low  power  densities  is  due  to  the  large 
decrease  in  pB op  caused  by  relative  increases  in  the  compressor  and 
fan  loads  Wlost  to  fuel  cell  power  output  Wgross  as  VVnet  decreases 
below  1500  W  as  illustrated  for  all  Tamb  in  Fig.  2. 


Results  from  the  system  model  identify  how  critical  parameters 
influence  both  overall  efficiency  and  water  balance  for  the  inte¬ 
grated  power  plant.  Baseline  conditions  for  user-specified  variables 
used  in  the  current  system  simulations  are  provided  in  Table  1.  For 
a  given  simulation  condition,  the  overall  system  efficiency  pth  is 
found  from  Eq.  (1 )  which  shows  pth  in  terms  of  fuel  cell  power  out, 
parasitic  loads,  Wlost,  and  enthalpy  input  of  the  hydrocarbon  fuel  in 

T^HC.in^HC.comb- 

VVnet  ncens  Acen  iVce ii  —  vi/jost  .  . 

77th  =  -r —  = - -r- - r -  (1) 

QhC  TTlFic.in  “HC.comb 

pth  can  be  broken  up  into  three  components:  the  combined  effi¬ 
ciency  of  the  fuel  processor  and  Pd-alloy-membrane  pw,  the 
efficiency  of  the  fuel  cell  itself  ppC,  and  the  efficiency  of  the  balance- 


Table  2 

Key  system  states  and  outputs  for  two  Wnet  baseline  conditions  from  Table  1 


Parameter 

Value 

Wnet  (W) 

2500 

5000 

Overall  system  efficiency,  ^th 

0.301 

0.268 

Individual  cell  voltage  (V) 

0.75 

0.70 

Gross  power  demand,  Wgr0ss  (W) 

3242.5 

6721.0 

Net  water  balance  (g  s_1 ) 

0.124 

0.071 

Reformer  outlet  temperature,  Tref,0ut  (K) 

952.7 

975.7 

WGS  membrane  reactor  conditions 

Inlet  reformate  temperature,  TWGs,ref,in  (K) 

822.2 

870.2 

Outlet  reformate  temperature,  TWGSirefi0ut  (K) 

859.8 

917.1 

Outlet  sweep-stream  temperature,  TWGs,Sweep,out  (K) 

889.9 

935.0 

Inlet  reformate  pressure,  PwGs.ref.in  (bar) 

7.31 

7.53 

Inlet  sweep-stream  pressure,  PwGs.sweep.in  (bar) 

1.31 

1.53 
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Fig.  3.  Parasitic  loads  for  LP  and  HP  compressors,  cooling  fans,  and  liquid  pumps  for 
the  liquid-fueled  PEM  fuel  cell  system  at  baseline  conditions  (see  Table  1)  at  three 
different  ambient  Tamb  and  Wnet. 


The  trends  in  rjth  as  a  function  of  Wnet  are  consistent  at  all  3 
Tamb,  but  also  as  shown  in  Fig.  2,  there  is  a  monotonic  decrease  in 
77th  with  increasing  Tamb-  This  is  due  to  an  increase  in  Wlost  due  to 
compressor  loads  rising  with  Tamb.  The  drop  in  ijth  is  approximately 
2%  points  for  each  20  °C  rise  in  temperature,  and  this  result  shows 
the  importance  of  stating  system  efficiency  values  at  a  given  Tamb- 

For  a  small-scale  fuel  cell  generator,  changes  in  total  parasitic 
loads  Wiost  (from  the  balance-of-plant  components  and  the  dc  to 
ac  inverter)  with  operating  conditions  determines  how  77Bop  influ¬ 
ences  the  overall  77th.  Fig.  3  shows  changes  with  Wne t  in  parasitic 
loads  at  three  Tamb  (for  the  otherwise  baseline  conditions).  The 
parasitic  loads  are  broken  out  into  groups  of  components— the  LP 
compressor  for  the  cathode  flow,  FIP  compressor  for  the  reformer  air 
supply,  cooling  fans  for  the  fuel  cell  radiator  and  exhaust  condenser, 
and  pumps  (including  all  liquid  pumps  and  the  anode  recircula¬ 
tion  compressor).  At  the  baseline  reformer  conditions  (S/C  =  1.6  and 
O/C  =  0.8),  the  LP  and  HP  compressors  are  the  chief  contributors  to 
the  parasitic  power  demand.  Both  compressor  demands  increase 
similarly  with  Wne t  and  Tamb-  The  total  power  for  the  fans  Wfan 
also  increases  with  Wne t.  Because  of  the  low  required  fan  speeds 
at  less  rigorous  operating  conditions,  the  Wfan  dependencies  only 
becomes  significant  for  Wne t  >  3000  W.  With  respect  to  Tamb,  the 
radiator  fan  power  demand  increases  with  Tamb  because  of  the 
larger  coolant  loop  heat  load  coming  from  the  fuel  cell.  However, 
the  power  demand  of  the  exhaust  condenser  fan  decreases  at  higher 
Tamb  because  of  the  smaller  heat  load  due  to  less  condensation. 
This  variation  however  does  not  significantly  impact  the  trends  in 
total  Wiost  due  to  the  fact  that  fan  power  at  these  conditions  is  still 
relatively  small. 

In  addition  to  77 th,  overall  system  water  balance  is  a  critical  mea¬ 
sure  of  the  system  performance.  Water  balance  is  measured  by  the 
difference  between  the  water  collected  from  the  exhaust  condenser 
and  the  anode  inlet  heat  exchanger  (condenser)  and  the  amount  of 
water  demanded  by  the  fuel  processor  and  the  sweep  steam  for  the 
low-pressure  side  of  the  Pd-alloy-membrane.  As  discussed  earlier 
and  illustrated  in  Fig.  1,  there  are  really  two  water  balances:  (1)  a 
clean  water  loop  involving  the  anode  inlet  condenser,  the  fuel  cell, 
and  the  sweep-steam  generator  and  (2)  a  fuel-processor-water  loop 
involving  the  autothermal  reformer  and  the  exhaust  condenser. 
Excess  from  the  clean  water  loop  can  be  supplied  to  the  fuel  proces¬ 
sor  water  loop,  and  thus  the  total  water  balance  will  be  presented  as 
a  single  value  in  these  results.  However,  because  the  anode  recir¬ 
culation  loop  is  essentially  closed  except  for  flow  across  the  fuel 


cell  Nafion  membrane,  water  recovery  from  this  loop  is  near  com¬ 
plete  and  independent  of  variations  in  Tamb.  Thus,  water  recovery 
from  the  burner  exhaust  condenser  is  most  critical  in  determining 
overall  water  balance  in  the  system. 

Fig.  2  also  shows  the  variation  of  net  system  water  balance  with 
Wnet  for  Tamb  of  10, 30,  and  50  °C.  As  is  consistent  with  other  system- 
level  studies  [10],  the  net  water  balance  was  most  sensitive  to  the 
outlet  temperature  of  the  exhaust  condenser  Texh-  For  the  current 
study,  Texh  is  set  to  Tamb  +  10°C  by  varying  the  fan  power  except 
for  conditions  (at  high  Wnet  and  lower  Tamb  which  leads  to  more 
condensation)  when  the  fan  flow  rate  reaches  its  limit  at  a  velocity 
of  3  ms-1  across  the  condenser  face.  For  the  range  of  conditions 
where  Texh  follows  Tamb  + 10  °C,  the  net  water  balance  follows  a  near 
linear  relationship  with  Wne t  as  illustrated  in  Fig.  2.  Such  linear 
behavior  implies  that  for  most  Tamb.  The  net  water  balance  grows 
either  more  negative  or  positive  with  Wne t  for  a  given  Tamb.  For  the 
cases  in  Fig.  2,  water  balance  remains  positive  for  Tamb  at  10  and 
30  °C  and  negative  for  50  °C.  At  high  Wnet,  the  fan  hits  its  flow  rate 
limits  and  Texh  rises  above  its  desired  valued  for  all  Tamb  <  50  °C.  The 
effect  of  this  fan  power  limit  on  the  water  balance  is  shown  in  Fig.  2 
where  the  net  water  balance  drops  below  its  linear  behavior  at  the 
high  Wnet  for  10  and  30  °C. 

At  a  given  Tamb*  the  overall  system  water  balance  magnitude 
increases  with  respect  to  Wnet  until  the  exhaust  condenser  air¬ 
flow  reaches  its  maximum  allowable  velocity.  Beyond  this  point, 
the  exhaust  condenser  cannot  recover  additional  water  as  the  sys¬ 
tem  demand  for  water  increases.  Thus,  further  increases  in  Wnet 
result  in  net  decrease  in  water  balance.  For  Tamb  equal  to  50  °C,  the 
current  system  does  not  meet  water  balance  at  any  load  conditions. 
In  this  case,  insufficient  water  recovery  in  the  exhaust  condenser  is 
due  to  the  high  Texh,  rather  than  the  air  velocity  constraint,  limit¬ 
ing  the  extent  to  which  the  burner  gas  cools  and  condenses.  Fig.  4 
illustrates  the  strong  dependence  of  overall  system  water  balance 
on  Tamb  for  half  and  full  power  conditions.  Fig.  4  also  indicates  the 
variation  in  r]th  with  Tamb  and  shows  a  drop  of  2-3%  points  when 
Tamb  rises  from  10  to  50  °C. 

For  the  current  system  configuration,  the  water  balance  is  posi¬ 
tive  for  Tamb  approaching  40  °C  at  full  power.  Near  40  °C,  the  water 
balance  shifts  from  positive  to  negative  water  balances  as  Wne t 
increases  due  to  the  exhaust  fan  reaching  its  flow  rate  limits  before 
Texh  cools  to  its  desired  value  of  Tamb  + 10  °C.  For  Tamb  >44°C,  the 
water  balance  becomes  negative  for  all  values  of  Wnet  for  the  given 
configuration.  The  difference  between  Texh  and  Tamb  can  be  reduced 
at  the  expense  of  condenser  size  and  fan  power  to  maintain  pos¬ 
itive  water  balances  at  higher  Tamb-  Also,  it  should  be  noted  that 


Fig.  4.  Overall  system  efficiency  and  total  net  water  balance  as  functions  of  ambient 
temperature  for  full-load  and  half-load  conditions. 
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Fig.  5.  (a)  Overall  system  efficiency  and  total  net  water  balance  as  functions  of  H2 
recovery  fraction  of  the  Pd-alloy-membrane  for  full-load  and  half-load  conditions: 
(a)  at  30  °C  ambient  and  (b)  at  50  °C  ambient. 

increases  in  the  reformer  O/C  ratio  relative  to  S/C  ratio  can  raise  the 
^amb  threshold  for  maintaining  positive  water  balance.  Finally,  fuel 
cell  Vcell ,  which  is  a  function  of  stack  size  and  polarization  curves, 
affect  water  balance,  since  the  non-linearity  in  the  water  balance 
with  respect  to  power  changes  is  due  mainly  to  the  non-linear 
change  in  the  stack  current  with  power  as  discussed  above. 

3.2.  Effects  of  fuel  processor  and  H2  purification  on  system 
performance 

For  this  fuel  cell  generator  with  a  liquid-fuel  reformer  and  a  Pd- 
alloy-membrane  purifier,  rjth  is  very  sensitive  to  both  the  fuel  cell 
operation  voltage  (through  p¥C)  and  the  membrane  H2  recovery 
fraction  (through  ^fp)-  The  membrane  recovery  fraction  is  defined 
by  the  fraction  of  total  FI2  that  crosses  over  to  the  low-pressure, 
purified  side  of  the  membrane.  As  p^p  decreases,  a  greater  total 
H2  flow  is  needed  from  the  fuel  processor  to  supply  the  same  FI2 
flow  to  the  fuel  cell  stack.  Increased  fuel  consumption  requires 
increased  air  and  water  flow  rates  to  the  fuel  processor,  along 
with  greater  parasitic  power  consumption.  Increases  in  membrane 
recovery  fraction  improve  pth  as  indicated  in  Fig.  5a  for  Tamb  =  30  °C 
and  Fig.  5b  for  Tamb  =  50  °C.  On  the  other  hand,  there  is  improvement 
in  the  water  balance  as  membrane  recovery  fraction  decreases.  This 
is  due  to  the  greater  availability  of  FI2  in  the  retentate  flow  and  thus 
H20  in  the  burner  exhaust. 


Fig.  6.  Critical  system  temperatures  as  functions  of  H2  recovery  fraction  of  the  Pd- 
alloy-membrane  for  VVnet  -  2.5  kW  at  ramb  =  30°C. 

The  recovery  fraction  of  the  Pd-alloy-membrane  not  only 
impacts  the  system  efficiency,  but  also  the  operating  temperatures 
of  many  of  the  system’s  components.  Fig.  6  illustrates  the  effect  of 
recovery  fraction  on  a  range  of  critical  high  T  locations  associated 
with  the  fuel  processor  and  the  exhaust  burner.  For  lower  mem¬ 
brane  recovery  fractions,  the  increased  H2  content  in  the  retentate 
flow  results  in  a  higher  outlet  temperature  from  the  exhaust  burner. 
Since  the  amount  of  heat  that  can  be  recovered  from  the  system  is 
limited,  any  extra  heat  must  be  rejected  to  the  ambient,  which  sig¬ 
nificantly  reduces  pth.  This  drop  in  pth  with  reduced  membrane 
recovery  fractions  is  exacerbated  by  higher  Tamb.  At  lower  Tamb  as 
shown  in  Fig.  5a,  sufficient  water  can  be  recovered  from  the  burner 
exhaust  condenser  to  achieve  a  net  positive  water  balance  across  a 
range  of  membrane  recovery  fractions.  By  contrast,  at  high  Tamb 
(50  °C)  as  illustrated  in  Fig.  5b,  the  higher  Texh  limits  the  water 
recovery  and  does  not  allow  for  improved  water  balances  with 
reduced  membrane  recovery  fractions.  Thus,  no  benefits  can  be 
achieved  by  reducing  membrane  recovery  fractions  at  high  Tamb 
operation.  The  trends  shown  in  Fig.  5a  and  b  indicate  that  the  design 
and  control  of  the  water-gas-shift/membrane  reactor  and  H2  puri¬ 
fier  is  central  to  achieving  an  efficient  and  operable  generator  when 
it  is  chosen  as  part  of  the  fuel  processing  sub-system. 

Additional  simulations  represented  in  Fig.  7a  and  b  investigate 
the  impact  of  varying  the  stoichiometry  of  the  liquid-fuel  reformer, 
while  keeping  the  sum  of  S/C  and  O/C  fixed  at  2.4.  This  fixed  sum  of 
S/C  and  O/C  is  prescribed  such  that  there  is  adequate  oxygen  con¬ 
tent  in  the  feed  to  avoid  coking  in  the  reformer  [23].  A  sum  of  S/C 
and  O/C  below  this  point  increases  the  risk  of  carbon  deposition 
on  the  reformer  catalyst,  thereby  degrading  performance.  Increas¬ 
ing  S/C  and  decreasing  O/C  significantly  improves  pth.  This  is  due, 
in  part,  to  the  greater  FI2  yield  from  steam  reforming  compared 
to  partial  oxidation  and  thus  a  rise  in  p^p.  However,  in  addition, 
reducing  O/C  also  reduces  the  flow  demand  from  the  HP  compres¬ 
sor  for  the  fuel  processor  air.  Reducing  the  high  power  demand  of 
the  HP  compressor  (relative  to  the  liquid  pump  for  the  reformer 
steam  supply)  increases  pBop  which  further  contributes  to  the  rise 
in  the  overall  pth.  Comparing  Fig.  7a  and  b  shows  that  the  influence 
of  high  Tamb  negatively  impacts  both  pth  and  overall  water  balance 
at  all  values  of  S/C  although  the  impacts  are  less  severe  at  half-load 
than  at  full-load. 

The  benefits  in  pth  with  increasing  S/C  are  not  offset  by  drops 
by  overall  system  water  balance  as  indicated  in  Fig.  7a  and  b, 
and  thus,  it  would  appear  that  S/C  should  be  increased  as  high 
as  possible  to  improve  system  performance.  However,  system 
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Table  3 

Breakdown  of  system  efficiency  at  full  and  half  power  for  a  range  of  S/C  and  7amb 


s/c 

Tamb  (K) 

VVnet  (kW) 

Off 

Ore 

??BOP 

0th 

1.4 

303.15 

2.5 

0.589 

0.598 

0.743 

0.262 

5.0 

0.589 

0.550 

0.707 

0.229 

323.15 

2.5 

0.589 

0.597 

0.725 

0.255 

5.0 

0.589 

0.542 

0.657 

0.210 

1.6 

303.15 

2.5 

0.649 

0.600 

0.772 

0.301 

5.0 

0.649 

0.555 

0.743 

0.268 

323.15 

2.5 

0.649 

0.599 

0.758 

0.295 

5.0 

0.649 

0.554 

0.733 

0.264 

1.8 

303.15 

2.5 

0.710 

0.601 

0.796 

0.339 

5.0 

0.710 

0.559 

0.772 

0.306 

323.15 

2.5 

0.710 

0.600 

0.786 

0.335 

5.0 

0.710 

0.558 

0.767 

0.304 

general,  the  increasing  endothermic  nature  of  the  reformer  with 
increasing  S/C  drops  the  temperatures  of  the  autothermal  reformer 
and  puts  a  greater  requirement  on  heat  recirculation  from  burner 
exhaust. 

For  autothermal  reformers  of  liquid  hydrocarbons  fuels,  studies 
have  shown  that  near  100%  conversion  of  the  hydrocarbons  and 
pyrolysis  intermediates  requires  reactor  temperatures  above  700  °C 
[24].  As  S/C  increases  much  above  the  baseline  value  of  1.6  (with 
O/C  decreasing  from  its  baseline  of  0.8),  the  existing  configuration 
cannot  sustain  such  high  temperatures.  While  clever  reactor  design 
with  sequential  injection  or  increased  heat  recirculation  may  be 
used  to  expand  the  range  of  S/C  available  for  the  system,  assessing 
the  feasibility  of  such  approaches  is  outside  the  scope  of  this  study. 


Fig.  7.  System  efficiency  and  total  water  balance  versus  fuel  reformer  S/C  (where 
fuel  reformer  O/C  is  set  to  2.4 -S/C)  for  full-load  and  half-load  conditions  (a)  at 
Tamb  =  30  °C  and  (b)  at  Gmh  =  50  °C. 

limitations,  including  the  challenge  and  expense  (in  terms  of  sys¬ 
tem  size  and  volume)  of  providing  low-cost  heat  to  the  reformer  as 
the  process  becomes  increasingly  endothermic,  make  it  beneficial 
to  provide  some  degree  of  oxygen  to  the  fuel  processor  to  increase 
the  exothermic  nature  of  the  fuel  reformer  reactions.  Fig.  8  shows 
how  variations  in  reformer  S/C,  with  corresponding  decreases  in 
O/C  impact  the  critical  high-temperature  locations  associated  with 
the  reformer,  WGS/membrane  reactor,  and  retentate  burner.  In 


Fig.  8.  Critical  system  temperatures  as  functions  of  S/C  (where  fuel  reformer  O/C  is 
set  to  2.4  -  S/C)  for  VVnet  =  2.5  kW  at  Tamb  =  30  °C. 


4.  Discussion 

Table  3  documents  values  of  ^pp,  ?7fo  and  r]  Bop  for  a  range  of  oper¬ 
ating  conditions.  For  the  range  of  operating  conditions  in  Table  3, 
trends  in  r]¥C  and  op  follow  each  other  with  changes  in  S/C,  Tam b, 
and  Wnet .  Since  decreases  in  r]¥C  result  in  greater  system  flow  rate 
requirements,  trends  in  ij¥C  are  mirrored  in  the  trends  of  r]B 0p  at 
these  operating  conditions.  However,  for  Wnet  below  1500  W,  the 
diminished  efficiency  of  the  LP  compressor  at  low  flow  rates  results 
in  an  opposing  trend  between  rj¥C  and  ?7Bop-  There  are  also  strong 
trade-offs  between  fuel  cell  stack  size  and  ?7FC,  and  between  stack 
size  and  r/BOp.  Therefore,  i)¥C  and  r/BOp  can  be  significantly  influ¬ 
enced  for  a  given  application  through  fuel  cell  stack  size,  which 
may  often  be  dictated  by  cost  targets  for  a  given  generator  applica¬ 
tion.  For  a  given  application,  the  preferred  stack  size  is  a  function  of 
the  Vcen  versus  i  curve  and  Wnet  demand  for  a  given  condition.  An 
increase  in  fuel  cell  size  not  only  improves  ij¥c  by  increasing  l/cell 
with  a  lower  i  but  also  increases  r/BOp  by  decreasing  VVlost  from  the 
balance-of-plant  due  to  the  lower  flow  rates  both  for  the  cathode 
and  fuel  processor  with  the  higher  Vcell  at  a  given  VVnet.  Thus,  for 
applications  where  cost  and  perhaps  system  size  is  not  as  critical 
as  ?7th,  large  fuel  cell  stacks  may  prove  to  be  a  viable  approach  for 
maximizing  system  efficiency. 

Identifying  other  areas  for  improving  overall  system  efficiency 
and  water  balance  is  one  of  the  underlying  goals  in  establishing  the 
simulation  tool  in  this  study.  It  is  clear  from  the  complex  behav¬ 
ior  of  the  system  water  balance  and  ijth  with  respect  to  Wnet  and 
to  a  lesser  extent  Tamb,  improving  system  performance  requires  an 
understanding  of  the  range  of  operating  conditions  and  load  pro¬ 
files  for  a  given  application.  Furthermore,  maintaining  as  possible 
a  positive  or  minimally  negative  instantaneous  water  balance  for  a 
stand-alone  generator  application  may  lead  to  changes  in  fuel  pro¬ 
cessor  stoichiometry  or  membrane  effectiveness  that  can  greatly 
influence  system  performance. 
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Clearly  from  Fig.  7,  one  of  the  key  ways  to  improve  r]th  primarily 
through  raising  ij¥C  is  increasing  fuel  processor  S/C,  but  this  can  only 
be  done  if  there  is  adequate  heat  recirculation  or  reactor  design  to 
ensure  near  equilibrium  fuel  reformer  outlets.  Another  means  of 
improving  performance  is  to  increase  ?7bop  by  lowering  parasitic 
loads.  From  looking  at  Fig.  3,  it  is  clear  that  the  two  compressors 
provide  the  largest  fraction  of  VVlost.  The  cathode  LP  compressor  par¬ 
asitic  load  can  be  reduced  by  running  a  lower-pressure  drop  stack, 
or  by  developing  compressors  with  improved  isentropic  efficien¬ 
cies  at  the  reduced  VVnet  conditions.  The  HP  compressor  load  may 
be  dramatically  reduced  by  finding  a  lower  pressure  H2  purification 
process.  This  may  be  achieved  by  ultra-thin  (<5  fxm)  Pd-alloy- 
membranes,  which  may  not  require  the  6.0  bar  differential  total 
pressure  across  the  membrane  and  thereby  reduce  the  contribu¬ 
tion  to  VVi0St  by  the  HP  compressor.  Another  means  of  reducing  the 
HP  compressor  load  is  to  use  a  preferential  CO  oxidation  (PROx) 
reactor  (as  proposed  in  other  studies  [8,10,25]  to  replace  a  Pd- 
alloy-membrane  for  H2  purification.  However,  PROx  reactors  add 
penalties  by  adding  additional  air  flow  requirements  (for  the  oxi¬ 
dation  stream)  and  by  reducing  the  fuel  cell  stack  performance  with 
reduced  H2  anode  concentrations  and  non-zero  CO  anode  concen¬ 
trations.  Because  of  the  concentrations  of  C02,  CO,  and  N2  in  the 
reformate,  PROx-based  systems  are  not  amenable  to  anode  recircu¬ 
lation.  While  such  systems  may  be  capable  of  utilizing  H2  exhausted 
from  the  anode  for  heat  generation  or  water  recovery,  they  can¬ 
not  utilize  it  in  the  fuel  cell  stack.  Studying  the  trade-offs  between 
lowering  the  pressure  of  the  fuel  reformer  and  reducing  the  fuel 
cell  stack  performance  with  a  PROx  reactor  is  left  for  a  further 
study. 

While  the  LP  and  HP  compressors  are  the  chief  parasitic  loads, 
the  cooling  system  also  demands  a  significant  power  at  high  VVnet* 
and  improvements  in  r]B op  may  be  achieved  through  optimized 
heat  exchanger  design.  Pressure  drops  across  the  exhaust  heat 
exchangers  which  determine  fan  parasitic  loads  are  highly  non¬ 
linear  (approximately  cubic)  with  respect  to  air  flow  rates.  This  in 
part  causes  the  increased  power  demand  of  both  cooling  fans  which 
contributes  to  the  consistent  drop  in  r/BOp  with  increasing  VVnet.  At 
high  Tamb  however,  parasitic  loads  for  system  cooling  are  some¬ 
what  offset  by  reduced  heat  transfer  requirements  in  the  exhaust 
condenser  due  to  lower  condensation  rates.  In  fact,  if  the  condenser 
approach  temperature  is  held  constant,  the  heat  removal  to  con¬ 
dense  large  amounts  of  water  in  the  exhaust  increases  with  lower 
^amb  (particularly  <20  °C).  This  increases  fan  power  for  the  con¬ 
denser  for  these  lower  Tamb  as  well.  This  provides  an  additional 
opportunity  for  higher  ^Bop  through  a  complex  fan  power  con¬ 
trol  algorithm  to  let  the  Texh  -  Tamb  difference  float  with  existing 
operating  conditions  in  order  to  optimize  fan  power  while  ensur¬ 
ing  positive  for  water  balance.  Resulting  improvements  in  r]Bop  and 
thus,  r]th  through  improved  fan  power  control  will  not  be  as  signif¬ 
icant  as  increases  from  reducing  the  compressor  loads. 


5.  Conclusion 

The  system  analysis  presented  in  this  study  can  effectively  be 
used  to  assess  the  performance  of  a  liquid-fuel-powered  PEM  fuel 
cell  system.  The  tool  provides  a  qualitative  basis  for  understanding 
the  sensitivity  of  overall  r]th  and  system  water  balance  to  the  differ¬ 
ent  component  parameters  in  the  fuel  cell  stack,  the  fuel  processor 
and  H2  purifier,  and  the  remaining  balance-of-plant.  The  system 
model  is  also  critical  for  identifying  operating  conditions  for  inter¬ 
nal  components  as  a  function  of  global  operation  parameters  such 
as  fuel  process  S/C  and  O/C  and  fuel  cell  stack  Vcell.  Thus,  the  model 
can  assess  which  approaches  for  improving  r]th  and  system  water 
balance  will  maintain  components  in  feasible  operating  tempera¬ 


tures  and  pressures.  The  quantitative  results  of  the  simulations  will 
thus  aid  in  designing  components  such  as  the  fuel  processor,  heat 
exchangers,  and  balance-of-plant  by  providing  the  designer  with 
the  range  of  possible  operating  conditions. 

The  current  study  indicates  that  stand-alone  operation  of  a  PEM 
fuel  cell  generator  with  a  fuel  processor  and  a  liquid  hydrocarbon 
fuel  source  is  feasible,  subject  to  constraints  on  ambient  condi¬ 
tions  and  operating  parameters.  For  the  system  detailed  in  the 
present  study,  the  overall  system  performance  is  highly  sensitive 
to  fuel  cell  stack  Vce\h  H2  purifier  recovery  fraction,  and  fuel  pro¬ 
cessor  S/C  and  O/C  ratios.  Generator  performance  also  varies  with 
the  environmental  variable  Tamb.  Isentropic  efficiencies,  particu¬ 
larly  of  compressors  for  air  supply  and  fans  for  cooling  along  with 
their  necessary  discharge  pressures,  dominate  the  losses  in  ijth  due 
to  balance-of-plant.  The  simulations  reveal  a  trade-off  at  higher 
Tamb  in  performance  between  r]th  and  net  water  balance.  For  the 
system  configuration  in  the  current  study,  water  balance  for  the  sys¬ 
tem  configuration  is  maintained  for  Tamb  up  to  35  °C  at  full  power 
demand  and  up  to  44  °C  at  50%  power  demand.  Maximizing  r]th 
favors  running  the  system  with  the  highest  fuel  processor  S/C  that 
still  allows  sufficient  internal  reformer  temperatures  to  sustain  high 
fuel  conversion  in  the  fuel  processor  and  the  exhaust  burner.  The 
revealed  design  trade-offs  show  the  importance  of  putting  even 
simple  component  models  into  their  system  context  to  understand 
how  they  influence  the  system  performance  as  well  as  how  the 
system  context  influences  their  design  and  performance. 
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